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. I ' Abstract 

5-1 . 

Olj' An ADM-like Hamiltonian approach is proposed for static spherically symmetric relativistic star con- 

figurations. For a given equation of state the entire information about the model can be encoded in a 
certain 2-dimensional minisuperspace geometry. We derive exact solutions which arise from symmetries 
^\. ' corresponding to linear and quadratic geodesic invariants in minisuperspace by exploiting the relation to 

JH , minisuperspace Killing tensors. A classification of exact solutions having the full number of integrations 

constants is given according to their minisuperspace symmetry properties. In particular it is shown that 
Schwarzschild's exterior solution and Buchdahl's n — 1 polytrope solution correspond to minisuperspaces 
with a Killing vector symmetry while Schwarzschild's interior solution, Whittaker's solution and Buch- 
dahl's n — 5 polytrope solution correspond to minisuperspaces with a second rank Killing tensor. New 
solutions filling in empty slots in this classification scheme are also given. One of these new solutions has 
a physically reasonable equation of state and is a generalization of Buchdahl's n = 1 polytrope model. 
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1 Introduction 

The subject of relativistic star configurations goes back to tlie very birth of general relativity itself when Schwarzschild 
gave his interior solution in 1916 |1|- Although a large body of literature now exists on the subject there still remain 
unsolved problems. For example, a theoretical understanding of some aspects of the interplay between the equation of 
state and the geometry of the model has only recently been achieved S . The stability of stellar models is another 
area in which only partial results are known concerning the relation between stability and the equation of state. 
In this paper we present a Hamiltonian formulation of the equations for spherically symmetric static equilibrium 
configurations. The Hamiltonian framework is then applied to show how known exact solutions are explained in 
terms of Killing tensor symmetries in a certain minisuperspace. A number of new exact solutions are also uncovered 
in this process. 

Spherically symmetric gravitational fields belong to a class of models in which the spacetime is foliated or sliced 
by a family of homogeneous hypersurfaces. The causal character of the slicing is timelike. These facts lead one to 
compare and contrast these models with another class of models, namely the spatially homogeneous models, also 
hypersurface homogeneous but with a spacelike slicing. From a mathematical point of view the difi'erence between 
star models and spatially homogeneous spacetimes is that the independent variable is a spacelike radial coordinate for 
stars while for spatially homogeneous models it is a timelike coordinate which is interpreted as a cosmic time function. 
The dramatic success and usefulness of Hamiltonian methods for exact solutions M, H , solution structure ||q] , quantum 
cosmology |7| and visualization of the minisuperspace dynamics |8| make it natural to ask if the Hamiltonian methods 
might also be useful for star models. In fact the ADM formalism works also for timelike slicings of the spacetime 
M . In place of the lapse function one must then introduce a radial gauge function. One consequence of this is that 
the Schwarzschild coordinates are unsuitable as a starting point for the Hamiltonian formulation. Instead they will 
correspond to a certain choice of momentum dependent radial gauge function. Also, as explained below one must be 
careful when introducing the matter terms to obtain a correct Hamiltonian. This paper is intended to serve two main 
purposes the first of which is to outline the Hamiltonian approach to star models in general while the second main 
purpose is to classify those exact solutions which are related to symmetries of the relativistic equations for static 
equilibrium. 

The importance of the Hamiltonian approach for exact solutions stems from the possibility to define a dynam- 
ical minisuperspace geometry which carries all information about the dynamics. We shall define such a dynamical 
geometry by using a special radial gauge function to be referred to as the Jacobi gauge. This is analogous to the 
Jacobi lapse which was introduced for spatially homogenous models in pCl, fll|. The resulting Jacobi geometry is 
conformally related to minisuperspace geometries corresponding to other choices of lapse function. Given such a dy- 
namical geometry the problem of finding and classifying exact solutions becomes a matter of studying geodesic orbits, 
a subject which has been treated from many points of view in the literature. There are two different mechanisms by 
which exact solutions can occur. A system of differential equations can be reduced either by means of a symmetry 
or by means of an invariant submanifold leading to symmetry solutions and submanifold solutions, corresponding 
to the classical notions of general and particular solutions respectively. Broadly speaking, symmetry solutions are 
more useful since they allow one to impose arbitrary initial conditions. Fortunately, they are also more amenable to 
systematic study by using the theory of Lie symmetries for differential equations. 

In the context of relativistic star models initial conditions correspond to matching conditions at the surface of 
the star and possibly inside the star to match the solution to one with a different equation of state above a certain 
pressure. For spherically symmetric models the matching condition is simply that the pressures of the two solutions 
must be equal at the boundary. In particular, the surface of the star is defined as the point of zeropressure and this 
automatically matches the solution to the Schwarzschild vacuum solution by Birkhoff's theorem M. In addition, a 
solution which goes all the way to the center of the star must be regular there. This is the well-known condition of 
elementary fiatness. It is discussed in sections in the context of the Hamiltonian framework. An alternative to using 
symmetry solutions is to assume a certain functional form for the radial dependence of a metric component (e.g. 
[hd) or for the energy density (e.g. |l3]). If the functional form involves some arbitrary parameters, then one may 
obtain physically reasonable submanifold solutions by adjusting those parameters in order to satisfy the boundary 
conditions. 

In the geometric formulation of the dynamicaj the Hamiltonian takes the purely kinetic form Hj = i J^^paPb 
where the Jacobi metric Jab is the metric of minisuperspace in the Jacobi gauge. In this picture the star configurations 
correspond to the timelike geodesies of the Jacobi geometry. The simplest symmetries of Hj are those which are 
related to Killing vectors of Jab ■ It is well-known that a Killing vector ^ gives rise to a linear constant of the motion 

■^The word dynamics is somewhat inappropriate here since we are dealing with static models. However, the radial variable 
will be seen to play a role which is identical from the mathematical point of view to the role played by the time in the closely 
related spatially homogeneous models. In particular the radial metric coordinate serves as the independent variable for the 
Hamiltonian. Lacking a more suitable term we will refer to the Hamiltonian equations of star configurations as "dynamics" . 



for Hj given by ^' pa- Similarly, constants of the motion may exist which are homogeneous polynomials of quadratic 
or higher degree in the momenta. In the quadratic case such constants of the motion have the form ^^^paPb where 
^ ^ are the components of a geometric object known as a second rank Killing tensor [M . In general, a constant of 
the motion of the n'th degree corresponds to a n'th rank Killing tensor. In particular a Killing vector is a Killing 
tensor of rank one. The precise relation between Killing tensors and Lie symmetries of the geodesic equations was 
given in lla]. 

In the ADM formulation the static star configurations are l+l-dimensional Hamiltonian systems. A criterion for 
the existence of Killing vectors and second rank Killing tensors for such systems was given in H . It was also shown in 
[0 how a Killing vector or a Killing tensor can be used to obtain variables which are adapted to the symmetry. Such 
symmetry adapted variables are needed for explicit integration of the equations of motion. There is a large number 
of exact solutions in the literature (see Wm for a partial list) some of them symmetry solutions and some submanifold 
solutions. One of our main objectives is to identify the symmetry solutions. 

We shall use the Killing tensor criterion of |5| to look for integrable relativistic star models. While a general 
solution to this problem is not available a present, even for Killing tensors up to second rank, we will be able to 
make substantial progress by using a certain ansatz for the symmetry adapted variables. The resulting family of 
solutions includes all (to this author's knowledge) known symmetry solutions as well as some new solutions. There 
is no guarantee that the symmetry solutions are physically reasonable. We shall therefore use the following rough 
criterion to separate solutions with a more reasonable equation of state from the obviously unphysical solutions. We 
consider the equation of state to be physical if the conditions p > 0, p > and < dp/dp < 1 are satisfied. The last 
of these conditions ensures that the velocity of sound is well-defined by Usound ~ i/dp/dp and that it does not exceed 
the speed of light. Note that an equation of state which is defined by some function p = /(p) may be physical at some 
densities but unphysical at others. Our condition is in fact rather strict and excludes for example the incompressible 
fluid used in Schwarzschild's interior solution. One of the new solutions found in this work is physical according to 
our criterion and generalizes Buchdahl's n— 1 polytrope solution fig ]. 

2 Hamiltonian formulation of static models 

We consider static spherically symmetric matter configurations where the matter is a perfect fluid described by some 
equation of state p = /(p). The stress-energy tensor is given by Tap = {p + p)uc,ui3 +pgai3 and the metric is required 
to satisfy the Einstein equations Gap ~ nTa/3 ■ The metric of spherically symmetric models is usually written as M 



where 



ds' = -e^^dt^ + e^^dr'^ + r^d^^ , (1) 

dn^ = d6'^ + sin^6'd(?I>^ , (2) 

is the metric of the 2-sphere. One reason why this form has been preferred in the past is that the Schwarzschild r 
variable is invariantly defined by its relation A = Aur^ to the area ^ of a 2-surface given by dt = dr = 0. One lesson 
to be learned, though, from the experience of spatially homogeneous models is that no preferred set of variables exists 
which is suitable for all situations. In particular, for the Hamiltonian method we need to modify the form (hi) of the 
metric. Using the analogy with spatially homogeneous models a natural starting point for a Hamiltonian formulation 
is to write the metric in the form 

ds'^ = -e^'^'^df^ + N^dR^ + e^^'df^^ , (3) 

where A'^ is a radial gauge function analogous to the lapse function in formulations with a spacelike slicing. For static 
models A is a function of the radial coordinate R. The relation to conventional variables is given be the relations 

P^ = iy , e"' = r , NdR = e^dr . (4) 

Therefore the conventional variables correspond to the specific choice of radial gauge, A^s, for which the radial 
variable is precisely the Schwarzschild r variable. We shall refer to this choice as the Schwarzschild radial gauge. In 
the terminology of [[17| this is an example of an intrinsic gauge choice. The square of the Schwarzschild gauge function, 
Ns'^, coincides with the "Schwarzschild correction factor" of Harrison et al.[[I8[ This factor determines whether or not 
there is a conical singularity at the center of the star model. It can be computed in an arbitrary radial gauge by the 
relation (a comma will denote differentiation throughout the paper) 

Ns=e^ = N{r,Ry' = Ne-^\l3^Ry' . (5) 



As for spacelike slicings the Misner variablesn /3 are defined by 

/?" = 1(2/3^ +/3^) , /3^=/3°+/3+, 

/3+ = i(/?^-/3^) ,/3^ = /30 - 2/3+ . 



(6) 



diagonalize the kinetic energy in the ADM Hamiltonian. It is also convenient to define the Taub radial gauge by 
N = Nt = 126^"° (l|].0 Using the Lorentz frame uj° = e^^'dt, lo^ = NrdR, lu^ = e^'de^, uj'^ = e^' sin6id0 the Einstein 
equations for the metric (H) become 

Goo = -^e-«^V - j,e-''\l3%^ + /?+ «) + e'^^^^^^^^ = .p , 

G,,=^ie-«^V-e-^(''''+'^")=Kp, (7) 

G22 = G33 = a^e-^-^V + ^e-^'3''(2/3"«« - /3+ ^) = Acp , 

where 

T = i [-(/?°«f + (/?,+«)'] . (8) 

Note that the above equations are given in the Taub radial gauge. Eliminating T and solving for the second derivatives 
yields 

The ADM Hamiltonian for a spherically symmetric static fluid should be given by bl [l7| 

H = -2N <^Vn"n'^(G<,^ - kT^h) , (10) 

where '^V = e^^ is the 3-volume element (as usual only deflned only up to a constant factor) and n" is the unit 
normal to the homogeneous hyp ersurf aces. It is constrained by the Einstein equations to the zero energy surface, 
H = 0. Using the prescription dlCj) the Hamiltonian takes the form 

i/= iAA(-po'+P+') + 24AAe^''"-^^%24KAAe«''%(/3°,/3+) , (11) 

where A/" is a relative radial gauge function deflned byn J\f — N/Nt and p is the pressure written as a function of 
the metric variables. As discussed above the Hamiltonian is constrained by the Einstein equations to the zero energy 
surface, H — 0. 

At this point it is not yet clear how the variation of the pressure functionp(/3'', /J"*") is to be done in the Hamiltonian 
framework. In fact it turns out that the pressure cannot be varied independently with respect to /j" and /9+. This 
can be understood if we write down the fleld equations obtained from the Hamiltonian (fill) and compare them with 
the Einstein equations. For that purpose it is convenient to use the corresponding Lagrangian in the Taub radial 
gauge JV — 1 

L = i [-W%f + W%r] - 246^^°-^"+ - 24.e«''" p(/3°, /3+) . (12) 

The equations of motion of this Lagrangian are 

dL d dL „„ 4,30-2/3+ 1-, 6/3" r,. 6/3" dp „ 

(13) 

For these equations to coincide with the Einstein equations (H) we see by inspection that the pressure function 
p{(3 ,/?"'') must satisfy the relations 



^We use capital indices (A, B, ...) for the Misner variables taking the values 0,+. (The third Misner variable f3~ is not 
needed in this context.) 

^This gauge was originally introduced by Taub for spatially homogeneous models. Misner called it the supertime gauge in 
that context. 

^ In the context of spacelike slicings Jantzen [E0|| uses the notation x = M~^ while Ashtekar [ pi| uses the term densitized 
lapse for N = 12A^. 



leading to 

dp dp , , 

This equation shows that p must be considered as a function oi u — fi^ — (5^ — 2/?"*" in order that Hamihon's 
equations as derived from the Hamiltonian ([0]) represent the Einstein equations correctly. We have thus shown that 
the Hamiltonian for spherically symmetric static models is given by 

H = \N{-p^ + P+') + 2^e^^°-^^^ + 24«AAe«'3" V{v) , (16) 

where we have introduced the notation Viv) — p for the pressure function. The momenta can be expressed in terms 
of the velocities by 

A number of remarks about the nature of the Hamiltonian (U&) are in order. Interestingly, it is the pressure rather 
than the energy density which appears in the matter term by contrast to the Hamiltonians of models with spacelike 
slicings such as minisuperspace models for example. Since the kinetic metric carries a Lorentzian signature we can 
characterize variables as being timelike, spacelike or null with respect to that metric. In particular /3^ — logr is a 
null variable which increases with distance from the center of the star. It follows that any timelike or null variable 
must also be increasing functions. Another useful feature of the Hamiltonian (hd) is the freedom to perform Lorentz 
transformations which preserve the kinetic energy. For some purposes we will also use null variables defined by 

i; = /3" - /3+ = i(/?' + 2/3^) , v = l3'' = i(-«; + 3v) . 

The pressure is related to the energy density by the conservation equation 

-^^P + P. (19) 

This is consistent with the variational constraint ( |l5[ ) on the pressure. Therefore the natural procedure to integrate 
the system is to first specify the equation of state p = p{p), then use the conservation equation to write p as a 
function of i/ and then insert that function into the Hamiltonian. The resulting Hamilton's equations should then be 
integrated with suitable boundary conditions. We shall require as usual that the metric is well defined at the center 
so that it satisfies the condition of elementary flatness there. In addition it is necessary that the pressure is zero at 
the surface of the star in order that the geometry can be joined smoothly to the Schwarzschild metric in the exterior. 
To see how to express the elementary flatness condition analytically in terms of our variables we write the spatial 
3-metric as 

TVs^dr^ + r^dn^ . (20) 

Elementary flatness at the center requires that the Schwarzschild radial gauge function A'^s (expressed in a general 
gauge in equation (Bl)) tends to unity llSl. Using null variables the condition becomes 

Ns-'=N-'e-^^l. (21) 

To analyze how the zero pressure limit appears in configuration space let ps be the value of p at the surface. It is 
useful to distinguish between the cases ps ^ (liquid-like equation of state) and ps = (gas- like equation of state). 
In the first case when the pressure tends to zero it follows from (hgt) that i/ must tend to a finite constant, Us. For a 
given equation of state and star mass m, let Qs{m) be the point in configuration space where p = 0. By varying m, 
the points Qs{m) form a curve in configuration space representing the star surface for the given equation of state. 
Since p is a function of the spacelike variable v it follows that the star surface curve is a timelike plane given by 
Us- In the case p^ = one cannot use (|l9|) to draw the same conclusions. 



The center of the star can be defined as the limit of zero 2- volume, e^'^ ^ or equivalently w = p'' + fJ 



— oo. 



Thus, this limit is lightlike in configuration space. Just as for minisuperspace models it is possible to perform a 
Penrose compactification scheme for the configuration space (cf. M). By the above remarks it then follows that the 
star center would be located at 9_ or i_. A more refined argument can be made to show that the center of the star 
is actually located at i_ but we will not pursue this point further in this paper. 
A global quantity of particular interest is the mass function m{R) given by 

r^ /■"' ,1 

m{R) = 47r / pR''^dR' = 47r / pe^" d/3'^ . (22) 



In particular the total mass of the star is the value of the mass function at the surface, M — m{Rs). The total mass 
can also be computed without integration by using the matching conditions at the surface for one of the two metric 
components which depend on the mass. To see how this works for the component e^'^ — e^"^ we use the notation 
Z = e^ and note that for the Schwarzschild exterior geometry we have Z = \ — IMjr. However, we cannot use this 
relation as it stands to calculate M because the value Zs, is not invariantly defined since we can always rescale the 
time variable by a constant factor. Taking that gauge freedom into account we have 



Z = B(\~ 2M/r) , 

where B is some arbitrary constant. To eliminate the gauge factor B we differentiate (t 
solve for M. Evaluating the resulting expression at the surface yields 



M = i 



1 + 



rZ r 



1 + 



r,R Z 
rZ.R 



(23) 
with respect to R and then 

(24) 



for Schwarzschild and arbitrary radial gauges respectively. For the second way of calculating the mass we make use 
of the metric coefficient A'^s^ = e^^ = (1 — 2M/r)~^ the last equality holding for the external geometry. The value of 
this coefficient is invariant since the Schwarzschild radial coordinate itself is invariantly defined. Solving for the mass 
gives 

M = i[r(l-iVs-')]^ . (25) 

Which of the expressions for the mass one uses is a matter of convenience. We note in passing that the two expressions 
together imply the identity 



(iVs')s = 1 + 



r.RZ 



(26) 



It is sometimes useful to adapt the system to the /3 variable. This is achieved by performing a boost with velocity 



1/2 

leading to the Hamiltonian 
H = 



/3" = -i=(2^" + ^^ 



Vz 



p-^ ^MP° + W^ 



Vs 



-Vsp' 



\Af{-po^+p+'') + 2'We 



2^313° 



+ 24AA^c^^'^^^^'3">(-^/3r 



(27) 



(28) 



Note that this mathematical structure is similar to that of a negative curvature Robertson- Walker model minimally 
coupled to a scalar field (see WA)- The only difference is in the values of the coefficients of /3 . We shall also need 
the Hamiltonian in the null variables (hsl) 



H = -2Npn,Pv + 24A/-e"'+'" + 24«;AAe'<"'+"'p(i/) 



(29) 



recalling that v is given by (Fly) . When integrating the solutions it is useful to write the Hamiltonian in terms of the 
velocities 

H = -\N-^w.RV,R + 24AAe"'+'" + 24KAAe"'™+"'p(//) . (30) 

Exact solutions can be either symmetry solutions with the full number of integration constants or submanifold 
solutions (also called special or particular solutions) which arises whenever an invariant submanifold can be expressed 
as an explicit functional relation on the phase space. Since the latter solutions axe only particular solutions they 
cannot in general be made to satisfy the physical boundary conditions. For the models we are considering in this 
paper the symmetry solutions are also general solutions since a single symmetry is always sufficient to integrate a 
2-dimensional Hamiltonian system. 

The analysis of the symmetries of finite-dimensional Hamiltonians constrained to a zero energy surface 



H{q,p)=T{q,p) + U{q)=Q 



(31) 



with a quadratic kinetic energy function T{q,p) — ^h^^ (q)pAPB is much facilitated by going to the Jacobi gauge 
J\f = J\fj = J2C/|~^ in which the problem is reduced to finding Killing vectors and Killing tensors in the Jacobi 
geometry 

dsj'' = 2|[/|/iAsd/dg'' . (32) 

The method has been described in a series of papers dealing with spatially homogeneous models (see W H, uM and 
references therein). 



The Jacobi metric for the static spherically symmetric fluid can be read ofi' from (ES) with the result (after 
rescaling by a numerical conformal factor) 

ds/ = -2Gdu;d^ , G = 6""+=*" + Ke^'-"'+'''>Viu) . (33) 

Any Killing tensor of this metric gives rise to a constant of the motion and a corresponding Lie symmetry of the 
Einstein equations. Whether or not such a symmetry exists depends on the form of the pressure function Viy). In 
the next section we will see how such functions can be found. 

3 Killing tensor conditions and pressure functions 

In this section we introduce an ansatz which makes it possible to determine equations of state for which the Jacobi 
metric (B3) admits a second rank Killing tensor. This will include Killing vector cases as well. Since we are using 
an ansatz we do not expect to find all possible solutions but it turns out, however, that we do recover a substantial 
fraction of the solutions which have been discussed in the literature for their physical or mathematical interest. 
Moreover, we show elsewhere fea] that the ansatz actually gives all possible Killing vector solutions . 

The Jacobi metric is (l+l)-dimensional and therefore admits three distinct types of Killing tensors as shown in 
[fe[|. The Killing tensor Kab may have either two non-null eigenvectors or a single null eigenvector. When there 
are non-null eigenvectors two subcases arise, one which corresponds to standard Hamilton- Jacobi separability and 
another one for which the equations of motion can be decoupled in suitable complex variables. The possibility of 
non-Hamilton-Jacobi separation is a consequence of the indefinite nature of the Jacobi metric. 

We now state the necessary and sufficient conditions for the existence of a second rank Killing tensor in the Jacobi 
geometry (p3|). A set of null coordinates (VK, V^) are adapted to a second rank Killing tensor of the 2-dimensional 
Lorentzian metric 

ds^ = -2G(W, V)AWAV , (34) 

precisely if one of the following conditions hold H] : 

G^ww = or Gvv = , null Killing tensor, linear decoupling, 

G^ww — G.vv = , non-null Killing tensor, Hamilton-Jacobi separation, 

G^ww + G^vv = , non-null Killing tensor, complex decoupling. (35) 

The Killing vector case is contained in the Hamilton-Jacobi separation case. As a first step we check if the original 
null coordinates {ui,v) defined in (llq), modulo a possible scaling, can be made symmetry adapted in this sense by 
some appropriate choice of pressure function Viy). To use the above conditions (p3) we define null variables {W, V) 
hy w = kW, V — V. This corresponds to a boost combined with a trivial scaling of the original Misner variables (3 . 
The metric function can then be written as 

Q ^ ^kW+3V ^ ^g3feW+3Vp(^) ^ (3g) 

Calculating the second derivatives and representing them in the variables {v, V), they can be displayed in the form 



e-^^^+^-'G.ww = fc'e-s^ + ^k^ [9P{u) - 3P'(z.) + ip"(i.)] 



U.VV 



^ 9e-6^ + 9^ [p(;,) + -p'l^u) + \T"{i^)\ 



It is apparent from these formulas that G^ww and G^vv cannot separately be set to zero. It follows that there is no 
null Killing tensor adapted to these coordinates. Complex decoupling is also impossible since the coefficient of the 
first term of the right hand side has the same sign in both expressions in ( p7| ) . Consequently that term cannot be 
canceled when taking the sum. We can, however, obtain Hamilton-Jacobi separation by taking the difference 

g-i2v+6.(^^^^ - G.vv) = {k" - 9)e-'*^ 
+«;[9(fc' - l)V{v) - 3{k^ + 3)V'(iy) + i(fc' - 9)r"{u)] . (38) 

It follows that the right hand side vanishes identically precisely when k^ — 9 and V {v) = 2P(i/). Comparison with 
equation (119) then leads to the equation of state p + 3p = 0. However, this is only a special case of the more general 
equation of state p -I- 3p = constant for which the general exact solution was found by Whittaker p3| and which is 
shown to correspond to a null Killing tensor in section W. 

Let us now consider more general transformations than the simple null scaling used above. To help us in the 
search for suitable conformal transformations we examine the structure of the matter potential for some equations of 
state of physical interest. For spatially homogeneous cosmological models with the usual equation of state p = (7— l)p 



Table 1: Expressions for pressure functions p = 'P{u) for some equations of state where po, po, P* and p* are constants. The 
corresponding fluid potentials are sums of exponentials for integer n. The expressions may be subjected to a gauge translation 
of the metric variable u = /3^. The type codes stand for polytrope-gamma law (PG), Tooper's equation of state (T), relativistic 
polytrope (RP), incompressible fluid-gamma law (IG), incompressible fluid (I), gamma law (G) and Buchdahl's n = 1 and n = 5 
polytropes (BPl and BP5). The velocity of sound for Buchdahl's n = 1 polytrope is less than the speed of light for pressures 
p < p*. Buchdahl's n = 5 polytrope satisfies p < p for p < p*. 

p = 'P{v) equation of state type 

{B + i)-"-i A"+i [e-(fl+i)''/("+l) _ l]"+i p = Ap'V("+i) + Bp PG 

{n + l)-"-lii'-"(e-'' - l)"+l p = {p/E')"/("+l) -l-np T 

(n + l)-ipo[e-("+l)"-l] p = po + np, [n = 1/(7-1)] IG 

PoCe"" - 1) P = PO I 

p^g-(n+i). p={7-l)p, [n = 1/(7-1)] G 

9p, (e^" - 1)2 p = 12,/p7p - 5p BPl 

6-'^7^p*e-''(l -e'')6 p = p^/^ /{7pl^^ - 6p^^^) BP5 

the minisuperspace potential is a sum of terms which are exponentials of linear combinations of the Misner variables. 
Integrating the conservation equation (119) one finds that this is also true for many of the equations of state which 
have been discussed in connection with relativistic star models. Most such equations of state are special cases of the 
relation 

p = Ap^/'^+i) + Bp , (39) 

which interpolates between a relativistic polytrope p = Kp^'^^'" of index n for small pressures and a gamma-law 
equation of state p = (7 — l)p for large pressures. The choice A = ^-"/("+''-) ^ B — n gives the relativistic ideal 
gas considered by Tooper VM. Other special cases of (139) are p = po + (7 — 1)~^P and an incompressible fluid 
p = Po- To this list can be added Buchdahl's n = 1 polytrope while Buchdahl's n — 5 polytrope is an example of an 
equation of state which is not of the form (By). The exponential representations of the above mentioned equations of 
state are given in Table pi Symmetry adapted variables for solvable spatially homogeneous models are either linear 
combinations of the Misner variables themselves or exponentials of null Misner variables. Also, as noted in ITJ the 
null variables (W, 1^) defined hy W — e"^ , V'"^^^ — e" are symmetry adapted for the Schwarzschild interior solution. 
In all cases the symmetry adapted potential turns out to be a polynomial in the symmetry adapted variables. These 
observations together lead us to the ansatz 

W = e-" = r , V'' = e" , (40) 

for the conformal transformation where a and /3 are constant parameters to be determined.n It will be shown that 

the above ansatz for the conformal transformation can be used to recover all (to this author's knowledge) known 

equations of state which lead to a general solution of the equations of motion. In addition the ansatz also leads some 

new exact solutions. 

In the null variables defined by the ansatz (WU) the Jacobi metric (BSh takes the form (upon rescaling by a constant 

factor) 

ds,7^ = -2GdWAV , 

(41) 

and we have defined Z — e^'^ = W~°'V''^'^ and f{Z) — p — 'P{\ log Z). Next we ask under what conditions on the 
parameters a and (3 and on the function f{Z) do the variables W and V become symm etry adapted with respect to a 
Killing vector or a second rank Killing tensor. For this purpose we use the conditions (pSl) and we therefore compute 
the second derivatives of the metric conformal factor from (Ull) with the result 

G,ww = (a - l)(a - 2)W'^-'W^^~^ + (3q - 1)(3q - 2)W'-'°'-'W''^-^ f(Z) 

-a(5a - 3)W''"- V^"- V(^) + a V"" V^"" V"(^) , 

G,vv = (3/3 - l)(3/3 - 2)I¥""V^'^"^ + (3/3 - l)(3/3 - 2)W^^"-V^'^"V(^) 

+W{W - 1)W^'"" V"- V(^) + W^W^-W"^-'^ f" {Z) . (42) 



* One could also use the slightly more general transformation obtained by replacing the first relation by the scaled form 
' = {kW)°' . However, it turns out that this additional freedom does not lead to any new solvable models. 



4 Null Killing tensor solutions 

In this section we determine the equations of state for which there exists a null Killing tensor in the minisuperspace 
Jacobi geometry. From the conditions (z)-(zzz) in section H we know that the Jacobi metric (hq) admits a null Killing 
tensor if one of G,ww or G vv vanishes identically. Our aim is to find all possible functions / for which this is the 
case subject to the ansatz (EtJ). To facilitate this analysis we replace the variable set (W, V) in m^ by {Z,V) with 
the result 

G,ww = i/i^/^-i-g^/.^o-.)/. ^(^ _ ^)(^ _ ^^y-60 
+ {3a ^ l)(3a - 2)Z-^f{Z) - a{5a - S)Z~\f'{Z) + a' f" {Z)] , 
G,vv = v^^^'^-^^'°'Z^^-'^^''^ [(3/3 - l)(3/3 - 2)1/-^'^ 
+ (3/3 -l)(3/3-2)Z-"/(Z)+ 9/3(3/3 -l)Z"V(^)+9/3^f"(^)] . (43) 

From the first of these equations we see by inspection that G,ww = precisely if either 

a = l, Z^f"{Z)^2Zf\Z)+2f{Z)=Q, (44) 

or 

a = 2, 2Z'^f"{Z)-7Zf'{Z) + l0f[Z)=Q. (45) 

In the first of these cases, ((mI), the two independent solutions are f{Z) = Z and f{Z) — Z^ leading to the pressure 
function 

p = -ae*'^^ + fce^"^ = -aZ'^ + bZ , (46) 

where a and b are integration constants. Using the conservation equation (119) the energy density is then given by 

p = 5aZ^ - 3bZ . (47) 

The equation of state itself can be written in the form 

a{5p + pf ^2b^{3p + p) . (48) 

It can be shown that p, p and dp/dp cannot all be positive for this equation of state. It is therefore unphysical 
according to our criterion. We conclude the examination of this case by noting that the special cases a — and 6 = 
correspond to the equations of state 3p + p — and 5p + p — respectively. The first of these equations of state was 
aheady encountered above as an exactly solvable case while in the second case the Jacobi geometry (^l|) is fiat. 

Moving on to the case (kq), two independent solutions are f{Z) = Z^ and f{Z) — Z^'^ . The pressure and energy 
density can then be expressed as 

p = -aZ^/^ + bZ^ , p = GaZ^''^ - 5bZ^ . (49) 

The equation of state can be written in the form 

a^{6p + pf ^b''{5p + p)\ (50) 

Again this an unphysical equation of state for which p, p and dp/dp cannot all be zero simultaneously. 
The second possibility for a null Killing tensor occurs if G,vv = which happens precisely if either 

13=1 , f"(Z) = , (51) 

or 

/3=|, 2Zf"{Z)+3f'{Z)=0. (52) 

In the first of these cases, (pll), we have p = —aZ + b and p = 3aZ — b leading to the equation of state 

3p + p = 2b . (53) 

The corresponding solution was found by Whittaker ^^. Recall that the special case fe = was shown in section 
y to admit a non-null Killing tensor. The corresponding Jacobi geometry therefore admits two inequivalent Killing 
tensors. Finally in the case given by (|5^) we recover Schwarzschild's interior solution with the equation of state 
parametrized as 

P-^a + bZ-^^\ (54) 

p = a . 
This concludes the discussion of possible null Killing tensor cases. 



5 Non-null Killing tensor solutions 



To investigate the existence of non-null Killing tensors we must take the sum and the difference of the second 
derivatives G^ww and G^vv and set the resulting expression equal to zero according to the conditions (|3q). After 
multiplication by a suitable non-vanishing function we obtain 

_^(l-a)/a^l2/3-3-3/3/<.(.^^^^ _ ^G,vv) = ~eB{Z)V'''^ + Z" ' "^ A{Z)V''^ 

-e(3/3 - l)(3/3 - 2)1/"^ -I- (a - l)(a - 2)^^/"^* = , (55) 

where we have introduced the notation 

A{Z) = af'[Z) - a(5Q - 3)Z-^f{Z) + (3a - l)(3a - 2)Z-^f{Z) , 
B{Z) = 9/?V"(^) + 9/3(3/? - l)Z-^f'{Z) + {3f3 - l)(3/3 - 2)Z-^f{Z) , (56) 

{ui,U2,u3,Ui) = {0,fi,-6f3,fi - 6/3) , 

with fi = 2 — 6/3 /a, e = 1 for Hamilton- Jacobi separation and e = — 1 for complex decoupling. To analyze when the 
above expression vanishes identically we consider the dependence on the variable V . If the exponents Vi {i — 1,2, 3, 4) 
are all different then we see that the form of the coefficient of V"* in ( p5| ) implies that either a = 1 or a = 2. In 
the first case we have V2 = 2(1 — 3/3) implying that /3 7^ 1/3. Inspection of the coefficient of V"^ then shows that we 
must have /3 = 2/3 implying 

(//i,i/2,i^3,i^4) = (0,-2,-4,-6) . (57) 

Equating the remaining two coefficients to zero we obtain 

A{Z) = f"{Z) - 2Z-'f'iZ) + 2Z-''f(Z) = , 
B(Z)=4f"(Z)+6Z~\f'(Z) = . 
However, it's easy to see that this pair of equations is incompatible. Considering now the second case q = 2, we have 

{1^1, !/2, 1^3, Ui) = (0, 2 - 3/3, -6/3, 2 - 9/3)) . (59) 

Comparison with the coefficient of V^^ shows that /3 — 1/3 leading to 

(//l,!/2,l'3,J^4) = (0,1,-2,-1) . (60) 

The conditions for the remaining coefficients then become 

A{Z) = 4f"{Z) - UZ~\f'{Z) + 20Z~\f(Z) = , 
B{Z) = f"(Z) = . 

This system also turns out to be incompatible. We have thus shown that at least two of the exponents Vi must be 
equal. 

Suppose first that vi = 1^2, i.e., a — 3/3. Then it follows that 

(ui, V2,vs, i^a) = (0, 0, -6/3, -6/3)) , (62) 

implying that there are now only two linearly independent functions of V in (|55[), namely V^ — 1 and V^^ = V~^^ . 
Setting the coefficient of the latter function equal to zero shows that one of the following two conditions must hold: 

a = l /3 = l/3, (I) 

Q = 2 /3 = 2/3 . (II) ^ ' 

As we shall see case (I) contains Buchdahl's n — 1 polytrope solution mM as a special case while case (II) corresponds 
to Buchdahl's n — 5 polytrope solution fed] as recently generalized by Simon M. The corresponding equations of 
state will be discussed in the next section. 

Going now to the remaining cases we note that ui 7^ 1/3 always holds since /3 7^ 0. Consider therefore ui — 1/4 i.e. 
fi = 6(3 leading to 

(//i,i.2,i'3,!^4) = (0,6/3,-6/3,0) . (64) 

There are now three linearly independent functions of V in (pa). Inspecting the coefficient of one of those functions, 
V^^ , we find that there are two possibilities, /3 = 1/3 or /3 = 2/3. The first of these values for /3 is not compatible with 
fi = 6/3. The second value implies a — —2. For the remaining two coefficients to vanish we then have the conditions 

-eB(Z) + 12Z'^ = -4e/"(Z) - 6eZ'^f'{Z) + 12Z'^ = , 
A{Z) = 4f"{Z) - 26Z-\f{Z) + 56Z-'^f{Z) = . 
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Again, this is an incompatible system. 

The next case is V2 ~ vs corresponding to /i = —6/3. We then have 

(i/i, !/2, 1^3, i^i) = (0, -6/3, -6/3, -12/3) . (66) 

There are again three linearly independent functions of V in (pSl). In particular, from the coefficient of V^* we find 
that either a = 1 or a = 2. The first of these values is incompatible with fi — —6/3 so we consider now the second. 
The remaining equations to be satisfied in that case are 

B{Z)=2[2f"iZ)+9Z-'f'iZ)+6Z-^f{Z)^0, 
ZA{Z) - 12e = 4Zf"{Z) - 14/' (Z) + 20Z'^f{Z) - 12e = , 

Again it can easily be confirmed that this system is incompatible. Note finally that we always have vs ^ v^. We have 
therefore exhausted the possibilities to satisfy (pa). 

6 Equations of state corresponding to non-null Killing tensors 

We now derive the equations of state corresponding to (I) and (II) of (163). The equation of state for case (I) will be 
discussed in some detail while we shall be content with just writing down the case (II) equation of state. In case (I) 
the condition ( p5| ) is reduced to 

- eB{Z) + Z^A{Z) = , (68) 

leading to the equation 

{Z^ - e)f"{Z) - 2Zf'{Z) + 2f(Z) = . (69) 

Two independent solutions of this equation are f{Z) — Z and f{Z) — Z^ + e. It follows that we can formally reobtain 
the equation state (kS) by letting e — > 0. However, as we shall see the physical properties of the equation of state is 
very much different when e is nonzero. We write the general solution as 

f{Z) = aZ'^ - 2bZ + at , (70) 

corresponding to the pressure function 

p = ae " — 2be " + ae . (71) 

It follows that pressure and energy density can be parametrized as 

p = aZ'^ - 2hZ + ae , 

(72) 
p = —baZ + GhZ — ae . 

The equation of state can be written explicitly in the form 

86^(2ae-3p-p) = a(4ae-5p-p)^ . (73) 

Setting a = we again have the equation of state 3p + p = 0. Setting 6 = gives the equation of state 5p + p — 4ae. 
It is convenient to write the equation of state as 



p = a{Z' - 2SZ + e) , 
p = a(-5Z^ + 65Z - ( 



2 . ..^ . (74) 



where S — b/a. The parameter a is then seen to be a scaling which does not matter for physical quantities such as 

dp /dp or the relativistic adiabatic index 7 = (1 + p / p){dp / dp) . Its absolute value is also irrelevant for inequalities 

such as p > and energy conditions such as p > 3p. However, its sign does matter and must be determined in 

each case. For example, suppose the equation of state is given by p = ax^ , p = ax corresponding to p — a^^ p^ . 

For a > this is a n = 1 polytrope while for a < it is clearly unphysical. Moreover, given any solution of the 

Einstein equations gap for any equation of state p = ^{x), p = '^{x), the scaling properties of the Einstein equations 

guarantee that the scaled metric a'^g^p together with the scaled equation of state p = a<^{x), p = a^{x) also solve 

the Einstein equations. In fact, this is the relativistic version of a mechanical similarity transformation. 

For the Hamilton- Jacobi case (e — 1) with 5 = 1 we recover Buchdahl's n — 1 polytrope solution ||lq]. The 

harmonic solution and the Hamilton- Jacobi solution with 5 7^ 1 have apparently not appeared in the literature 

before. We now take a first look at the physial properties of these equations of state. The velocity of sound is given 

by the formula 

An A _ 7 

(75) 



2 dp 5- Z 



dp ^Z - 3(5 
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Since Z > hy definition the condition < dp/dp < 1 can only be satisfied if 5 > and 

^<Z<S. (76) 

The energy condition 

p + p^4aZ{-Z + 5)>0 , (77) 

then shows that we must have a > 0. Considering now the condition p > we first note that if e — —1, then that 
condition is nowhere satisfied in the region given by (|76|). Thus the harmonic solution is nowhere physical according 
to the adopted criterion and we shall not consider it further. 

For the Hamilton- Jacobi solution, e = 1, we find that ii S < 1 then p > automatically while if 5 > 1 we have 
the further condition 

Z <5~ ^/S^~l . (78) 

Compatibility of ( |78[ ) with (|7q) now leads to a restriction on 5 given by S/3 > \/S^ — 1. Since 5 > this is equivalent 
to 

^ < 3^ ~ 1.06 . (79) 

Summarizing our analysis so far, the conditions p > and < dp/dp > 1 are both satisfied if e = 1 and the conditions 

0<5<-^ , f<Z< fiS) , (80) 

are satisfied and where f{S) is defined by 

Examining now the condition p > 0, we find that the polynomial —5Z^ + 6SZ — 1 must have two distinct real roots 
(excluding the Z — constant case). It follows that 5 > \/5/3 ~ 0.745. Combining the condition p > with ( pO| ) gives 
the final result that the equation of state is physical according to our criterion in the range 

Y < ^ < 9iS) , (82) 

where 



„,„.|,3.WM^,/=,„.s;, „3, 



The requirement 25/3 < g{5) further restricts the parameter values to the range 

! < <5 < ^ . (84) 

Furthermore the limit p ^ can only be reached if 5 > 1. Therefore, if 5 < 1, this equation of state cannot be used 
all the way out to the star surface. One way comparing the equation of state with S j^ 1 with Buchdahl's original 
equation of state is to compute the ratio Ptj p, between the pressure and energy density in the limit Wsound ~^ 1- The 
result is 

J. = 165^^ ■ ('') 

This shows that the ratio is smaller for (5 > 1 and bigger for 5 < 1. A larger ratio can be interpreted as corresponding 
to a stiffer equation of state. Another way of examining the physical nature of the equation of state is to calculate 
the relativistic adiabatic index 

7= -T- ■ (86) 

p dp 

It can be shown that 7 increases with 5. This is also an indication that equations of state with larger 5 are stiffer or 
less compressible. 

We turn now to case (II) of (|63[). The pressure function is given by 

p = / (Z) = aFi (Z) + 6F2 {Z) , (87) 

where 

Fi(Z) = Z' + f eZ+1 , 

1/2/ ,73 I 1 £^5^2 



F2{Z) =Z'''''{Z'' + 15eZ^ + 15Z + e) 
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Using the variable Y = Z^'^ — e" the equation of state for the Hamilton- Jacobi case (e = 1) can be written in 
parametric form as (following Simon [H]) 

P = {6Y)-' [-p+(l + Yf + p-(l - Yf] , 

p = p+(l + yr + p-(l-y)^ 

where 

p+ = — = (a + 6fe) , , ^ 

' (90) 

p- — —-^(a — bo) . 

The subcase p+ = is Buchdahl's polytrope of index n — 5 p5| . The general case p+ 7^ was recently discussed by 
Simon M. 

The equation of state for the harmonic case, e — —1, can be written as 

p = (6y)-i Ke[po(l + jy)'] , P = 9m[po(l + i^)"] , (91) 

where po is a complex constant. We have not found any region where this equation of state is physical. 

7 Solving the configuration equations 

In this section we outline the procedure for integrating the Hamiltonian configuration equations for equations of 
state which correspond to a Jacobi geometry which admits a second rank Killing tensor. We begin by expressing the 
Hamiltonian (BOl) in terms of the symmetry adapted variables defined by the ansatz (H) 

H = ~lW,RV,R + 24al3[W"-W'^-^ +KW''"-W^^-^f{Z)] , (92) 

where the relative radial gauge is given by 

AA = al3{WV)-^ = 12q/3W'"+3'')/(3'')z-i/'^''^ , (93) 

while the absolute radial gauge becomes 

N = NrAf = 12q/3H^(='"~'')/ V^"""'''' . (94) 

We are interested in how the elementary fiatness condition (gl|) restricts the integration constants as we perform 
the integration. For that purpose we need to express the Schwarzschild radial gauge (H) in the symmetry adapted 
variables 

^ J_/3-lp^"(l-3/3)/(3/3)^(2-3^)/(6/3)^^ ^ ( ' 

In the symmetry adapted variables the Hamiltonian (|92[) has the form 



H = -^W.rV.r + A{W)V + B{W). (96) 

The procedure to obtain the general solution for this type of Hamiltonian has been outlined in [h7|. The fact that 
the potential in (Bq) is linear in one of the null variables, (V), leads to a decoupled equation of motion for the other 
null variable, (W), 

W,RR = 2A{W) . (97) 

This equation has the first integral 

'2\ 



E=\ {W,R^) - 2 / A{W)AW . (98) 

By integrating this equation we can express W as a function of R. Inserting this into the Hamiltonian constraint 
results in the linear first order equation for V given by 

V,R = {2/W,r)[A{W)V + B{W)] . (99) 

The solution of this equation is 

V = 2W,fl / B {W^rY^ AR . (100) 



Rather than deriving the metrics for all the possible models we shall give a few examples which illustrate the 
technique. First of all we show how to calculate the metrics for Schwarzschild's exterior and interior solutions in 
the Hamiltonian framework. After that we calculate the metric for the generalization of Buchdahl's n = 1 polytrope 
solution. In the final subsection we discuss how the remaining models can be integrated. The equations of state 
which correspond to minisuperspace Killing tensors are summarized in table H. 
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Table 2: Perfect fluids corresponding to Killing tensor solutions. In the first two columns a "-" means that the value of a or /3 
is arbitrary. The parameters a, b, p± are real constants while po is a complex constant. The notations Y = e'^ and Y± = 1 it y 
are also used. For the case a = 2, /3 = 2/3, it is not possible to write down an equation of state in closed form. Instead 
the expressions for p{Y) are given in the fourth column. The notation "H-J" in the last column stands for a Hamilton-Jacobi 
Killing tensor type. The Hamilton-Jacobi case with (a, /?) = (1, 1/3) reduces to a Killing vector solution when a = b. 



a 


P 


p = p{Y) 






equation of state 

[p = p{Y) in the last two rows] 


Killing tensor type 


1 


- 


-oY* + bY"^ 






a(5p + p)2 = 262(3p + p) 


null 


2 


- 


-aY^ + bY'^ 






a4(6p + p)5 =65(5p + p)4 


null 


- 


1/3 


-aY'^ + b 






3p + p = 2b 


null 


- 


2/3 


-a + bY-^ 






p = a 


null 


1 


1/3 


oY-^ - 2bY^ 4 


- a 




(5p + p - 4a) 2 86^ 
2a - 3p - p a 


H-J 


1 


1/3 


ay4 _ 2by2 - 


- a 




{5p + p + 4a) 2 8fe2 
2a + 3p + p a 


harmonic 


2 


2/3 


{6y)-i (-P+ 


rt + p- 


ye) 


p+Yf+p.Y^ 


H-J 


2 


2/3 


(6y)-i SRe [pc 


,{l + iYf 


Sm[po(l + iy)5] 


harmonic 



7.1 Schwarzschild's exterior solution 

To illustrate how the Haniiltonian method works for the static models we first derive the Schwarzschild vacuum 
solution. In this case there is no matter term and choosing a = 1, /3 = 1/3, the Hamiltonian (M) becomes 



H = -^W,rV,r + 8 



(101) 



The problem is now trivially solvable. In fact the Jacobi geometry is obviously flat in this case. The Hamiltonian 
constraint gives W,r — C and V^r = 16/C where C is a separation constant. Since r = e'^ = W it follows that R — 
C^V (up to a gauge translation) and V = 16C~^{R- Ro). The radial gauge is given by iV = NtJV = AW^^^V'^^'. 
Using the relations to the metric variables, e — W V , N — 16WV , one then arrives at the usual form of the 
Schwarzschild metric with the mass given by M = CRq/2. 



7.2 Schwarzschild's interior solution, p = Ps 

The equation of state is parametrized according to (B4J) as p = f{Z) — —a + bZ~^'^ and p = Ps 
13 = 2/3 and we set a = 1 implying W = r and Z — W~^V^. Referring to (M) we also have 

A{W) = 16(1 - KaW'^) , B{W) = 16KfeW^/^ . 

Applying the procedure outlined above we have the following first integral 

32(W~ iKaW"^) . 



a. Further we have 



(102) 



E = ^W,R^ 



(103) 



The next step is to use the condition of elementary flatness (pih to determine the value of E. It follows form (Bq) 
that W^R can be expressed in terms of the Schwarzschild radial gauge as 



W,R = SW'-^^Ns 



(104) 



Elementary flatness requires that TVs -^ 1 as r = W — * . Th is implies that W^r = at the center of the star. It 
then follows that E = for regular models. The equation (|l03|) can then be integrated further in terms of elementary 
functions. Alternatively one can use an intrinsic variable ]17|] such as VK or a function of W as coordinate in place of 
R. In fact, it is the intrinsice variable choice W — r which gives us the usual form of the metric. In practice, what 
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we must do to use an intrinsic variable is to change the integration variable in (100). Using W as the new integration 
variable that equation reads 

V = 2W,R / B{W) (Wm)-^ dW . (105) 

Performing the integration gives 

V = ialV^/^ (36-^ - kkJi - ^KaW^ j , (106) 

where K is an integration constant. Using Z = W~^V^ gives one of the metric coefficients 

Z ^ e^" = lb" kK Ud - ,Jl^^\^^A , (107) 



where D = (naK) ^. Finally, using (104) to compute the remaining metric coefficient yields 

e"^ = iVs' = 64W(W^,fl)-' = (l-iK/9,r')"' . (108) 

Using (psl) we also obt ain t he mass from this expression as M = Kp^Ts^ ■ The physical meaning of the constant D 
comes from evaluating (|l07|) at th e sur face. This gives D = (1 — irepsra^)^'^ where we have used Zs — (b/a)^. These 



relations together with (|107D and (108) give the standard form of Schwarzschild's interior metric as given in ||l4|. The 



constant b corresponds to a gauge scaling of the time coordinate and consequently has no physical meaning. 

7.3 Buchdahl's polytrope of index one and its generalization 

In this section we calculate the metric components for the solution corresponding to case (I) of (p3). This case 
generalizes Buchdahl's n = 1 polytropic solution to a fluid which is no longer gas-like in the limit p = but rather 
has a non-zero density there, ps 7^ 0. Since a — 1 and P = 1/3 we have Z — V/W and the symmetry adapted null 
variables are W — e^ and V = e"^". Specializing the Hamiltonian (M) to this case yields 

H^-^W,RV,R + »[l + Ka{V^ -25WV + W'^)\ . (109) 

Separation is achieved by introducing the symmetry adapted non-null variables (T, X) by 

T=\{W + V)^\e''\l + e'f'') = \r{l + Z) , 
X = \{W^V) = \e^\l^e^^''') = \r{l~Z) . 
In terms of these variables the Hamiltonian takes the manifestly separable form 

H^\{-T,R^ +X,R^)+S[l-2tia[5-l)T'^ + 2tia[5 + l)X^] , (111) 

leading to the separated equations 



(110) 



(112) 



Tm + 32Ka(5 - l)r^ =K+IG , 

X^R^ + 32Ka{S + 1)X'^ = K , 

where K is the separation constant. It is obvious from the second of these equations that K > 0. It is also clear 
that the solutions T(R) will have d iffer ent functional forms for the three cases 5 < 1, 5 = 1 and 5 > 1. Since T is a 



cyclic variable in the Hamiltonian (111) when 5 = 1 it follows that the original Buchdahl n = 1 polytrope solution 
corresponds to the existence of a Killing vector in the minisuperspace Jacobi geometry. 

As explained in section M, when 5 < 1 the equation of state has no zero pressure limit. Although such an equation 
of state could still be useful for high pressures we restrict the following discussion to the cases with S > 1. The 
solutions are given by 

(4{coshOiR-R-) , for<5 = l , 

^ ' \4{coshOLO--^sm[Lu^{R-R-)] , for (5 > 1 , (113) 

X{R) =4(sinhC)a;+"'sin[u;+(i?-i?+)] , 



where (^ — arcsinh(VA/4), lj± — •y'32aK((5 ± 1), and the parameters R± are int egra tion constants. Let us now 



consider the possible range of values which T and X can assume. It is obvious from (110) that T > 0. Also, it follows 



from (p2h that Z < 1 throughout the interior of the star implying that X > 0. It is also clear from (IIC) that both 



T and X vanish at the center. Next we use the translational freedom in the radial coordinate to make it vanish at 
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the center, Re — 0. The condition r(0) = then implies R- — for 5=1 and that uj-R- is a multiple of n when 
S > 1. Further X{0) — implies that uj+R+ must be a multiple of n. Since a/3 = 1/3 > it is clear from (M) that 
the radial gauge function is positive and hence R increases outward from the center. To ensure the positivity of T 
and X when _R > it follows that iiJ+R+ and uo-R- must both be even multiples of vr. But then we can just as well 
set R± = without loss of generality. 

We now consider the elementary flatness condition (Ell) at the center of the stellar model. In this case the 
condition can be transformed to 

lim iWflZ'/^ = 1 . (114) 

For the above solution we have {T^r)c = 4cosh^, (X,_r)c = 4sinh(^ leading to {W^r)c = 4e'' and {V^r)c = 4e^''. The 
central value of Z is given by the indeterminate expression Zc — {V/W)c which by the previous argument is seen to 
be given by Zc = e~^^. Using these central values we find that the elementary flatness condition is already satisfied 
without restrictions on the remaining integration constant (. Restricting attention now to the case S > 1, the final 
form of the solution is therefore given by 



T{R) = 4{coshQ{u}-y'^ sm{iO-R) 



(115) 



Referring to (fm) and ( [Tq ) we must require for consistency that 25/3 < Zc < Zb leading to a restriction on the 
separation parameter 

25/3 < e-'^^ < 5 - y/P~i . (116) 

The central values of the pressure and energy density can be inferred from ([7^ as functions of S and (^. Using ( |94| ) 
to calculate the radial gauge we find that the metric for this model is given by 

ds^ = -Zdt'^ + mZ-^dR^ + W^dn'^ , (117) 



where Z = V/W = (T - X)/{T + X), W = T + X aad T and X are given by (yj). The physical properties of this 
model will be given elsewhere. 

7.4 Other models 

We begin this subsection by discussing how the remaining null Killing tensor models can be integrated. Whittaker's 
solution with equation of state given by (|53| ) has a structure which is similar to Schwarzschild's interior solution. The 
Hamiltonian is again linear in V and thus the solution can be integrated in an analogous manner. In this case the 
condition for elementary flatness does not restrict the integration constant E to be zero but the resulting equations 
can nevertheless be integrated in terms of elementary functions in the Schwarzschild radial gauge. This leads back 
to the form for the solution given by Whittaker |23| . For the null Killing tensor cases corresponding to the equations 
of state (Uq) and (BOl) the roles of W and V are reversed. This gives a Hamiltonian which is linear in W leading to a 
decoupled equation in V . Although the integration proceeds in an analogous way it is less straightforward to impose 
the elementary flatness condition. This is because that condition is expressed in terms W . However, an argument 
can be made using the Hamiltonian constraint to express the flatness condition in terms of V. These solutions cannot 
in general be expressed in Schwarzschild coordinates. 

The second Hamilton-Jacobi solution corresponding to the equation of state (p9|) ca n b e integrated in much the 



same way as Buchdahl's n — 1 polytrope and its generalization were treated in section 7.3, Buchdahl's equation of 
state corresponding to p+ = tends to a Newtonian n = 5 polytrope in the zero pressure limit. In fact, like the 
Newtonian n — 5 polytrope, the star has an inifinite radius. However, the more general Simon solutions do not seem 
to have this property and a closer examination of the physical properties of this class of solutions would be of interest. 
Although the harmonic solutions found in this work do not seem to of physical interest we briefly indicate the 
method by which they can be integrated. The decoupling variables in these cases are complex and given by Q = W+iV 
and its complex conjugate Q = W — iV . Using these variables together with the pressure functions ( |74| ) (with e = 1) 
and (M) the Hamiltonian (B3) becomes explicitly separated into a Q-dependent and a Q-dependent part respectively. 
The procedure to obtain the metric variables then parallels that for the Hamilton-Jacobi case. 

8 Concluding remarks 

We have shown that a number of known exact solutions for static star configurations can be described in a unified 
setting using an ADM- like minisuperspace geometrization of the dynamics. New solutions have also been found in this 
process illustrating the power of the Hamiltonian formalism. The remark by Kramer et al. in 1980 (|14|, p. 130-131) 
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that "A Hamiltonian formulation ... is not well-adapted to searching for exact solutions" reflected a prejudice of the 
time which should be dispelled by now. Indeed, the traditional view towards exact solutions is that, as stated by 
Schutz in 1985 when discussing Buchdahl's n — 1 polytropic solution, that "Finding such [exact] solutions is an art 
which requires the successful combination of useful coordinates, simple geometry, good intuition, and in most cases 
luck" (Pql, p. 263). In the light of the present work as well as earlier work on spatially homogeneous models (see 
m/n and references therein), the opposite view seems perhaps more appropriate. At least those exact solutions which 
correspond to symmetries and also many submanifold solutions [E7| can certainly be derived by systematic methods. 
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